Self-consistent study of Anderson localization in the Anderson-Hubbard model in two 

and three dimensions 



Peter Henseler and Johann Kroha 
Physikalisches Institut, Universitdt Bonn, Nussallee 12, 53115 Bonn, Germany 



oo 

o 

O 

(N 

O 

Q 

oo 



i 

C 

O 

o 



> 

cn 

06 
o 
oo 
o 



X 



Boris Shapiro 

Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 

We consider the change in electron localization due to the presence of electron-electron repulsion in 
the Anderson-Hubbard model. Taking into account local Mott-Hubbard physics and static screening 
of the disorder potential, the system is mapped onto an effective single-particle Anderson model, 
which is studied within the self-consistent theory of electron localization. We find rich nonmonotonic 
behavior of the localization length £ in two-dimensional systems, including an interaction-induced 
exponential enhancement of £ for small and intermediate disorders although £ remains finite. In three 
dimensions we identify for half filling a Mott-Hubbard-assisted Anderson localized phase existing 
between the metallic and the Mott-Hubbard-gapped phases. For small U there is re-entrant behavior 
from the Anderson localized phase to the metallic phase. 

PACS numbers: 71.10.Fd, 73.20.Fz, 72.15.Rn, 71.30.+h 



I. INTRODUCTION 



The research on interacting, disordered systems is 
one of the central topics in today's condensed-matter 
physics. In particular, the experimental signatures of a 
two-dimensional (2D) metal-insulator transition in dilute 
disordered electron systems^ absent in the single-particle 
theory of Anderson localization,—^ have triggered many 
theoretical research activities. As a prototype model 
for the interplay of strong correlations and random- 
ness, the Anderson-Hubbard model has been intensely 
studied^ - — While in electron systems the bare Coulomb 
interaction is long ranged, inducing the additional ques- 
tion of the range of the effective interaction in a disor- 
dered system, the Anderson-Hubbard model assumes an 
on-site interaction which may grasp many of the salient 
features of the electron-electron interaction. Disordered 
many-body systems with tailor-cut short-range interac- 
tions have most recently been realized experimentally as 
cold atomic gases in random potentials ^ Despite these 
efforts, even in systems with short-range interaction the 
existence of a metallic ground state in d < 2 dimensions 
has remained elusive. Recent numerical works showed 
that the presence of interactions can at least suppress 
significantly the localizing effect of the disorder 4"—^ 

One of the ideas proposed to explain the delocalizing 
effect is screening of the random potential by the inter- 
action, discussed controversially in the literature. A de- 
tailed analysis of this screening effect is the subject of 
this paper. For that purpose, our approach starts from 
the atomic limit of the Anderson-Hubbard model. In the 
presence of a Hubbard interaction U, the on-site energy 
of a particle depends on whether this site is occupied by 
another particle or not, leading to an interaction-induced 
renormalization of the distribution of on-site energy lev- 
els. Below we argue that this effective disorder distribu- 
tion of the atomic limit will still provide a good descrip- 



tion of the static interaction effects when a finite hop- 
ping amplitude is included in the full Anderson-Hubbard 
model as long as the fluctuations of the on-site energy lev- 
els are large compared to the kinetic band energy (hop- 
ping amplitude). In particular, this will hold for arbitrar- 
ily large Hubbard interaction even when the localization 
length becomes large. As seen below, the latter may oc- 
cur even for relatively strong disorder near the Anderson 
transition in three dimensions and as an interaction ef- 
fect in two dimensions. By this reasoning, the Anderson- 
Hubbard model is reduced to an effective single-particle 
Anderson disorder model^ with renormalized level distri- 
bution, for which we calculate ensemble-averaged single- 
particle and transport properties at temperature T = 0. 
Clearly, inelastic as well as virtual interaction effects are 
neglected by this approach. For the present purpose 
they may, however, be unimportant due to the vanish- 
ing quasiparticlc relaxation rate^ at the Fermi energy 
in 2D and three-dimensional (3D) disordered systems for 
T = 0. 

In a previous workji 7 - we already presented an analyti- 
cal study of this approach, exploiting an exact relation^ 
between the localization length and the ensemble- 
averaged single-particle density of states (DOS) in one di- 
mension. In accordance with numerical results obtained 
for the Anderson-Hubbard model^ - — we could demon- 
strate that weak interaction reduces the effective disorder 
(screening) while a strong interaction effectively enhances 
the localization, corresponding to a hopping suppression 
by interaction (Mott-Hubbard physics). 

In this paper, we extend our analysis by applying 
the self-consistent theory of Anderson localization^ - — 
which allows for a quantitative analysis of the effective 
single-particle Anderson disorder model in a broad pa- 
rameter regime, particularly in two and three dimensions. 
Despite the simplifications made in our approach, we find 
good agreement with recent numerical studies, especially 
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in the Anderson localized regime of one and two dimen- 
sions. Our analytical approach therefore allows for a crit- 
ical assessment of some of the conclusions drawn from the 
numerical results. Our results indicate that screening of 
the disorder seems to be the most relevant physical mech- 
anism for the interaction-induced delocalization effect in 
the disorder localized regime of the Anderson-Hubbard 
model. 



II. ATOMIC-LIMIT APPROXIMATION 

We consider the Anderson-Hubbard model for fcrmions 
at zero temperature on a hypercubic lattice in d dimen- 
sions with lattice spacing a. It is described by the Hamil- 
tonian 



H = H + H 



kin 



= E i £ i - h c *W - * E c V c j° 



<i,j>,a 



^E'v 



(1) 



c\ a (c ia ) are creation (destruction) operators of a fermion 
at site i with spin <r and n ia — c\ a c ia . t is the 
nearest-neighbor hopping amplitude, U > is the 
on-site repulsion, and \i is the chemical potential. The 
on-site energies {s i } arc assumed to be independent 
random variables with a box probability distribution 
p(e) = 0(A/2 — |e|)/A, where the disorder strength is 
parametrized by the width A. The lattice filling, i.e., 
the average particle number per lattice site (summed 
over spin), will be denoted by p. 

In the atomic limit, t = 0, the energy of a particle 
on site i depends on the occupation number rij of that 
site and can be obtained in the paramagnetic phase by 
shifting the bare on-site energy according toii 
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U 



if fx 

if e, > n 



TT ( each with 

U <*i<» (prob.ofl ) 



The probability distribution pa(s) of these renormalized 
energy levels is, hence, identical to the (averaged) spec- 
tral density of the Hubbard two-pole Green's function 
in the atomic limit (shifted by the chemical potential). 
Examples of this renormalized distribution are given in 
Fig. [T]and in Ref . fl7l There it was shown that for U <C A 
the variance of Pa(s) as compared to p(e) is reduced. 
On the other hand, it is evident from Eq. (|2|) that for 
U > A the support of the distribution Pa(s) splits into 
two disconnected intervals, leading eventually to a Mott- 
Hubbard gap in the averaged DOS of the model [Eq. |0J)] 
below. 

When a finite hopping amplitude t > is switched on, 
the particles become dclocalizcd from a single site (with 



finite or infinite localization length £). The hopping in- 
duces quantum fluctuations of the occupation numbers 
m, and the on-site energies are renormalized by self- 
energy corrections of leading relative order 0[(t/A) 2 ]. 
However, since the average occupation number on each 
site is essentially determined by the minimization of the 
local electrostatic energy, the atomic limit approxima- 
tion will still capture the essential static physics of the 
Anderson-Hubbard model (disorder screening and Mott- 
Hubbard physics), as long as t <C A, for arbitrary U 2& 
In particular, this remains valid even for arbitrary local- 
ization length £ (Ref. l23h since the charge density may 
vary on the scale of a lattice spacing, independently of the 
size of £. With these assumptions the Anderson-Hubbard 
model is mapped onto an effective single-particle Ander- 
son disorder model. 
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(3) 



with the renormalized on-site energy distribution pa(s). 

In d = 1, 2 dimensions, as well as in d = 3 dimensions 
for sufficiently strong disorder, all particles described by 
Hamiltonian © arc exponentially localized^ The decay 
of their wave functions ip(r), in the limit r — > oo, is gov- 
erned by the localization length £. To analyze the effect 
of the repulsive interaction on the localization, we study 
the U dependence of £. For a first qualitative estimate 
we used in Ref. Qj] the relation^ (from now on we choose 
units where a = t = 1) 

oo 

er 1 = J N{e)\og\E-e\ de 

— OO 
OO 

« J PA (e) log \E - e + m| de, (4) 

— oo 

valid in one dimension, where £i is the wave-function 
decay length, N(e) denotes the disorder-averaged DOS, 
and E is the particle energy measured relatively to the 
chemical potential. The second approximate equality in 
Eq. (|3]) holds for A 3> t. For sufficiently large disor- 
der, Eq. ((3]) becomes also a good approximation in d > 1 
dimensions^ In Ref. [l7| we found that for all lattice 
fillings p the localization length of a particle at the Fermi 
level, £i, is a nonmonotonic function of U with a pro- 
nounced maximum at intermediate U. 

III. SELF-CONSISTENT TRANSPORT THEORY 

In order to extend our analysis of the effective single- 
particle system [Eq. l[5))] to two- and three-dimensional 
systems as well as to parameter regimes with large local- 
ization lengths and to the calculation of general trans- 
port properties, one must go beyond the restrictions of 
one dimension and strong disorder implied by Eq. ((¥]). 
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Therefore, we study the system within the self-consistent 
theory of Anderson localization^ - — This theory consti- 
tutes a resummation of the most divergent (Cooperon) 
contributions to the irreducible particle-hole vertex, lead- 
ing to a self-consistent equation for the dynamical dif- 
fusion coefficient. The theory was originally developed 
for weak disorder and was extended to arbitrary disor- 
der later on i 21 i 25 i 26 By comparison with direct numerical 
diagonalization results^ it was demonstrated that for 
the non-interacting Anderson model this theory yields 
quantitatively correct results for the phase diagram of 
Anderson localization in d = 3 and for the localization 
length in d = 1, 2, 3 dimensions (with exception of the 
critical regime). 

For the interacting case with short-range, instanta- 
neous interaction U, it is believed on general grounds 
that Fermi-liquid theory remains valid in the presence of 
disorder. Then the diffusion pole structure of the density 
propagator is preserved for particles at the Fermi energy 
(E = 0) for temperature T = 0. This was recently shown 
with the use of Ward identities^ 3 - to hold strictly at least 
when disorder- and interaction-induced self-energy con- 
tributions may be taken to be additive. 

In the formulation of Refs. Hfl and[26| the self-consistent 
equation for the diffusion coefficient D(lo,E) reads 



D{u,E) = D {E) 



2ImE( J B) 



d d k 



[ImG(E)] 2 D (E) J {2n) d 



d d k' ^lmG k (E){ImG k ,(E)} 2 



(2tt) c 



(k + k') 2 -iuj/D{uj,E) 



(5) 



where Dq(E) is the bare diffusion constant 
I f d d k 



D (E) 



ImG(E) J (2ir) d 



(v k -qr[ImG k (£)n6) 



and q is the unit vector in the direction of the transport. 
The disorder-averaged retarded single-particle propaga- 
tors are given in terms of the self-energy as 



G k (E) = [E + fi-^-ZiE)}- 1 , 



(7) 



and G{E) = J ^G k (E). e k = -2£ti cos(fc 4 ) is the 
dispersion and v^- = Vk£k is the group velocity. 

In Eq. ([5]) the diffusion pole structure of the integral 
kernel holds strictly for Q = |k + k'| < 2n/£, where I 
is the elastic mean-free path. For larger Q the integral 
kernel does not vanish but behaves in a non-singular way. 
Therefore, the momentum integrals in Eq. §5§ must not 
be cut off for Q > 2ir/£ but extend over the complete first 
Brillouin zone of the lattice . 25 ! 26 Furthermore, one has 
to keep in mind that not only the single-particle Green's 
function but also the particle-hole propagator obey lat- 
tice periodicity with respect to their momentum argu- 
ments. In particular, the particle-hole propagator, which 
enters into the integral kernel of Eq. Q ) 20 ' 26 is lattice 
periodic with respect to the center-of-mass momentum 



of the particle-hole pair. However, in Eq. ([5]) this peri- 
odicity is not explicit since the diffusion pole form of the 
kernel arises from a hydrodynamic expansion for small Q 
and lo. To restore the lattice periodicity in the transport 
properties, the subscript p in the denominator of the ker- 
nel of Eq. ([S]) implies a shift by a reciprocal-lattice vector 
so as to keep the momentum argument k + k' within the 
first Brillouin zone. 

The localization length £ for particles at the Fermi en- 
ergy is defined as the exponential decay length of the 
density correlation function in the static limit , 20 i 28 



lim 



-iuj 



(8) 



where in the localized phase the diffusion coefficient is 
purely imaginary to first order in lo and vanishes linearly 
for lo — > 0. 

For the evaluation of the theory, we first cal- 
culate the disorder-averaged single-particle quantities, 
i.e., the self-energy and the local Green's func- 

tion G(E), within the well-known coherent-potential 
approximation 2 ^ (CPA). The CPA is known to interpo- 
late disorder-averaged single-particle quantities correctly 
between the limits of weak and strong disorders, neglect- 
ing only exponentially rare disorder configurations (Lif- 
shitz tails of the DOS), and provides quantitatively re- 
liable results for disorder-averaged single-particle quan- 
tities over the complete parameter ranged 6 - It is defined 
in connection with Eq. (J7]) and the renormalized level 
distribution pa(z) by the self-consistent relation, 



depA(e) 



e - Z(E) 



l-[e- X(E)]G(E) 



0. 



(9) 



After the single-particle quantities are determined, the 
diffusion coefficient D(lo — > 0, 0) or the localization 
length £, respectively, is calculated by solving numeri- 
cally Eqs. © and © with Eq. (jSJ. The two static in- 
teraction effects, disorder screening and Mott-Hubbard 
gap formation, are incorporated in the single-particle 
transport theory through the quantities, Gu_(E = 0) and 
H(E = 0), determined by the renormalized distribution 
Pa{s). 

In our static treatment of the interaction term, the 
hopping of a particle at the Fermi level was assumed to 
happen on a background of immobile particles. In the 
case of singly occupied sites, the spin of these particles 
was considered to be randomly distributed. Therefore, 
in the absence of the site-energy disorder, i.e., A — ► 0, 
Eq. © reduces to the Hubbard III approximatio n 30 ! 31 
and, consequently, can be understood as an average over 
both spin and site-energy disorder. Implications and re- 
strictions on the applicability of our approach will be 
mentioned in the discussion of our results in the follow- 
ing sections. 

In d > 1, the evaluation of the 2D integrals in Eq. ([5]) 
is numerically costly. Taking advantage of the period- 
icity of the integrands, this problem can be overcome 
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FIG. 1: CPA density of states (solid lines) and renormalized 
site-energy distribution pa{e) (dotted lines) in d — 2 at half 
filling for A/t = 8 and (a) U/t = 0, (b) 4, (c) 8.5, and (d) 12. 
Energies are measured relative to the Fermi level. 



by rewriting the factors in the integrand as Fourier se- 
ries and using the convolution theorem. While for weak 
disorder the evaluation of the Fourier series is numer- 
ically still somewhat costly because of the pronounced 
peak structure of the Green's functions, the Fourier se- 
ries converges quickly for larger disorder. The solution of 
the self-consistent equation [Eq. (5])] is easily performed 
on a single desktop computer. 



IV. RESULTS AND DISCUSSION 

The transport theory described above allows one to 
analyze in detail the static disorder screening effect of the 
Hubbard repulsion U in the Anderson-Hubbard model in 
arbitrary dimension. 



A. Two-dimensional systems 




FIG. 2: (a) Localization length £ at the Fermi level (E — 0) 
in d — 2 as a function of U for A/t = 8 and two differ- 
ent band fillings. The arrows indicate the positions of the 
maximum of £ as obtained from Eq. Q for a ID system 
(Ref. [HI). At the Mott transition for p = 1 [vanishing JV(0) 
due to Mott-Hubbard gap formation], £(U) vanishes expo- 
nentially as ln£(?7) ~ — 1/2JV(0) because of the logarithmic 
divergence of the Cooperon integral in Eq. (J5j in o? = 2. The 
inset shows the critical interaction strength, U C (A), for the 
Mott transition at half filling, (b) Localization length £ at 
the Fermi level as a function of U in d = 1 and 2 at half 
filling, p — 1, for A/t — 4. The exponential enhancement of 
&d(U) ~ exp[l/A^ ff ] in d = 2 due to interactions is clearly 
seen while in d — 1 the dependence is comparatively weak 
(Refs.dJandUl), £ iD ~ [A eH (U)]- 1/2 (see text). 



Figure Q] shows the 2D DOS N(E), computed within 
CPA from Eq. ©, for a fixed disorder strength A and 
different values of the repulsion U. The figure clearly 
exhibits the regime of interaction-induced screening of 
disorder for small values of U, characterized by a nar- 
rowing of the disorder-averaged DOS with increasing U 
[Figs. Hlja) andQJb)]. It also shows the regime of hopping 
suppression for large U, where the Mott-Hubbard gap 
gradually develops with increasing U [Figs. [llb)-[ljd)]. 
The crossover between the two regimes occurs roughly 
at U « A. Note that the present static approximation 
does not describe the Kondo-like quasiparticle resonance 
at the Fermi level that would be induced by dynamical 
processes in high dimensions near half filling^ There- 
fore, it does not capture the typical many-body effects at 
the metal-insulator transition and inside the correlated 



metallic phase, known to be important when A — > 0. 

Figure [2] shows the generic behavior of the localiza- 
tion length at the Fermi level, £(U), as a function of 
the repulsion U for (a) strong and (b) intermediate dis- 
orders A. The most salient feature seen in Fig. [2] is 
the nonmonotonic behavior of the localization length 
with a pronounced maximum at an intermediate value 
< < A. Within the one-dimensional (ID) or strong 
disorder™ approximation [Eq. |QJ] , can be calculated 

as U 6 pa C/ ? (1D) = A[y/l + 3p(2-p) - l]/3^l As seen in 
Fig. [21 this provides even in d = 2 an excellent quantita- 
tive estimate for the results of the self-consistent theory 
not only for large disorder, as expected, but also for inter- 
mediate disorder. The nonmonotonic behavior is repro- 
duced by numerical methods for finite-size systems, i.e., 



by quantum Monte Carlo (QMC) simulations^ and sta- 
tistical dynamical mean field theory (DMFT)^ with the 

maximum of £(£/) occurring almost precisely at U^ 10 ^ . In 
the case of the QMC results^ the nonmonotonicity can 
be inferred from the dependence of the finite-temperature 
conductivity on U. In Ref. @ the on-site energies were 
calculated as poles of the atomic limit Green's func- 
tion of the Anderson-Hubbard Hamiltonian, resulting in 
precisely the same Hamiltonian as our effective model 
[Eqs. §2§ and J3J)]. In Ref. d this Hamiltonian was then di- 
agonalized numerically exactly for finite-size systems ac- 
cording to the statistical DMFT approach, and the local- 
ization length was extracted from the disorder-averaged 
inverse participation ratio (IPR) using the definition, 



£ := (IPR) 



-l/d 



(10) 



We find good agreement with the results of Ref. H for 
all parameter values available up to a factor of order 
unity. This factor might be attributed to the slight dif- 
ference in their and our definitions of £ [Eqs. |(5J) and 
(fTU)) ] . respectively. This agreement lends additional sup- 
port to the quantitative correctness of the results of the 
self-consistent transport theory within the atomic limit 
approximation. 

In our semianalytic theory the nonmonotonic behav- 
ior of £ is easily traced back to the two competing 
interaction-induced effects: screening of the random po- 
tential and DOS suppression due to a Mott-Hubbard gap. 
Both effects are already incorporated in the effective dis- 
tribution pa(s), as seen in Fig. [1] under an increase in 
U . the effective disorder is initially reduced, leading to 
an increase in £. For large U, however, the formation of 
a Mott-Hubbard gap implies a reduction in the DOS at 
the Fermi level and a broadening of the disorder distribu- 
tion and, hence, a reduction in £ with increasing U [Fig. 
Ha)]. 

In Ref. the interaction strength was fixed and the 
IPR as a function of the disorder strength was studied. 
That work was based on a local unrestricted Hartree- 
Fock treatment of the interaction and subsequent ex- 
act diagonalization of the resulting effective disordered 
single-particle Hamiltonian. Our results show a non- 
monotonic behavior of the IPR with a pronounced sup- 
pression (i.e., increase in £) for intermediate disorder, as 
seen in Fig. [3f a) , in good agreement with the finite-size 
data of Ref. Our theory indicates that this results 
from an intricate competition of Anderson localization 
and correlation- induced disorder screening effects: for 
small increasing A the IPR increases from zero due to 
Anderson localization in d = 2 for arbitrarily weak dis- 
order. As A increases further to intermediate values, the 
effective disorder potential, derived from Eq. @, is re- 
duced by the screening effect until A«[/, resulting in a 
suppression of the IPR. Finally, when A > U, the on-site 
repulsion U cannot induce a further significant change in 
local occupation numbers as A grows; hence the effective 
disorder potential and the IPR increase. 
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FIG. 3: (a) Inverse participation ratio £~ 2 as a function of 
the disorder strength A in d — 2 for U/t = 4 at half filling 
(p = 1). The exponential suppression of the IPR for small 
disorder and the absence of the nonmonotonic behavior in the 
non-interacting case (U — 0) are also shown (dashed line), (b) 
Magnitude of the derealization effect as a function of disorder 
strength in d — 1 and 2. 



Because of the logarithmic divergence of the Cooperon 
integral in Eq. ([5]) for d = 2, the disorder screening ef- 
fect can induce an exponentially large effect as seen in 
Fig.OJb). To demonstrate this screening- induced dereal- 
ization effect quantitatively, we show in Fig. EJb) the ra- 
tio £(t^)/£({7 = 0) of the maximal interaction-enhanced 
localization length and its non-interacting value. The 
effect is especially strong for weak disorder in = 2, 
where £ depends exponentially on the screening-reduced 
disorder strength A e g of the effective Anderson model 
[Eqs. © and ©] £ ~ exp[l/A 2 eS ]M 

In Refs. iand0 also the observation of a metallic phase 
was reported. The existence of such a phase is not possi- 
ble in an effective single-particle model^ and would con- 
sequently be beyond our approach. However, the good 
agreement of our results with the finite-size data of these 
works and the observation of the exponentially enlarged 
localization length, which exceeds the largest system size 
used in the numerics, suggests strongly that the infinite- 
size extrapolations erroneously indicated a true metallic 
state in d = 2^ 
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FIG. 4: (Color online) 3D phase diagram in the (U, A) plane 
for different lattice fillings p. The metallic phase (M) is to the 
lower left while the insulating phase (I) is to the upper right 
of the solid mobility edge curves for p — 1. 00, 0.75, and 0.50, 
respectively. The dashed curve shows the critical interaction 
strength for the Mott transition at half filling. The dotted line 
is an extrapolation of the mobility edge for p = 1 toward small 
disorder, where the evaluation of the self-consistent Eq. © 
becomes numerically costly. The corresponding Anderson, 
Mott-Hubbard, and Mott-Hubbard assisted Anderson insulat- 
ing regions at half filling are marked by AI, MHI, and AI^ M ', 
respectively. 



The remarkable quantitative agreement of our theory 
with the results of numerical calculations wherever com- 
parison is possible suggests that the essential physics of 
localization in the 2D Anderson-Hubbard model is cap- 
tured by the assumptions of static disorder screening and 
Mott-Hubbard gap formation. 



For the half-filled case (p = 1) a further increase in U 
eventually leads to a suppression of the averaged DOS 
at the Fermi level, N(0), because of the gradual forma- 
tion of a Mott-Hubbard gap (dashed line) . Consequently 
the system undergoes first a transition from the metal- 
lic phase to a (Mott-Hubbard assisted) Anderson insula- 
tor with a reduced but finite A^(0) just before the DOS 
vanishes at the dashed line and the Mott-Hubbard insu- 
lating phase is entered. This intermediate phase is not 
likely to be seen in DMFT studies of the proble m 12 i 13 
because within DMFT the formation of a Kondo reso- 
nance and the concatenated unitarity sum rule for N(0) 
in the metallic state prevent a Mott-Hubbard-induced 
suppression of N(0) and lead to a complete screening of 
the disorder potential . 12 ' 13 It remains to be seen whether 
in the physical 3D disordered systems Kondo physics or 
pscudogap formation dominates. 

Away from half filling, a pure Mott-Hubbard transition 
is not possible because iV(0) remains non-zero. There- 
fore, for weak disorder the metallic phase extends out 
to arbitrarily large U. However, for large U (U/t> 10) 
increasing disorder A does induce an Anderson metal- 
insulator transition, and this transition occurs at a sub- 
stantially smaller value of A than for small and interme- 
diate U/t < 7 because of a Mott-Hubbard reduced DOS 
at the Fermi level (pseudogap). It is interesting to note 
that this interplay of Anderson localization and Mott- 
Hubbard physics away from half filling leads to a step-like 
behavior of the phase boundary, as seen in Fig. 2J con- 
necting smoothly to the phase boundary of the half-filled 



CONCLUSION 



B. Three-dimensional systems 

In d > 2, the non-interacting Anderson model de- 
scribes a disorder induced metal-insulator transition, 
where extended and localized states are separated by the 
mobility edge4 Hence, the question about a possible dis- 
order reduction by interaction gets even more relevant. 
We have calculated the phase diagram of localization in 
d = 3 using the atomic limit approximation [Eqs. @ and 
([3])] with the self-consistent transport theory of Anderson 
localization [Eqs. ©-©; sec Fig.g]. 

The figure shows, for different values of the lattice 
filling p, the disorder screening effect: above the crit- 
ical disorder value for Anderson localization without 
interaction 3 ^ A c /t w 11.7, the increase in interaction 
leads to a rc-cntrance from the Anderson insulating into 
the metallic phase, i.e., the metallic phase is extended 
to larger values of A. This behavior was also observed 
in a recent DMFT study^ of the Anderson-Hubbard 
model where a site-dependent self-energy correction was 
used. Contrarily, generalized DMFT studies based on 
site-independent averaged self-energy corrections do not 
describe the screening effect^ 



We have presented a self-consistent study of the static 
disorder screening effect and of Mott-Hubbard gap for- 
mation, induced by the local repulsion in the Anderson- 
Hubbard model. Both of these effects are represented 
by an interaction-induced rcnormalization of the effec- 
tive distribution of on-site energy levels and a subsequent 
mapping of the Anderson-Hubbard model onto a single- 
particle Anderson model. While this mapping is exact 
in the atomic limit, we have argued that it still provides 
a good description of the static screening effect for fi- 
nite hopping amplitude even when the localization length 
£ is large. The localization properties of the effective 
single-particle Anderson model were then treated by the 
self-consistent theory of Anderson localization. 19 ! 26 We 
found rich behavior of the localization length in two di- 
mensions and of the phase diagram in three dimensions 
due to an intricate interplay of disorder screening and 
Mott-Hubbard physics in the different regions of param- 
eter space. 

Despite the technical simplicity of our approach, it 
yields good agreement with numerical studies^ - — of the 
same problem for two-dimensional finite-size systems, in- 
cluding the non-monotonic dependence of £ on both the 
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interaction strength U and the disorder A. At the same 
time, we found in d = 2 an exponential interaction- 
induced enhancement of £ for weak and intermediate 
values of A although a true metallic state is not possi- 
ble within our effective single-particle theory. The good 
agreement of our results with the numerical finite-size 
calculations, and simultaneously, our prediction of a large 
but finite localization length show that the indications of 
a true metallic state, found in some works by infinite-size 
extrapolations of the numerical data, are not conclusive. 
These indications might rather be due to the fact that 
the largest system sizes were still smaller than the local- 
ization length of the infinite system. 

For three-dimensional systems we found that for weak 
Hubbard interaction U the disorder screening effect re- 



sults in a re-entrant behavior from the insulating to the 
metallic phase while for large U , disorder and Hubbard 
interaction cooperate to form a Mott-Hubbard assisted 
Anderson insulating phase (with finite density of states 
at the Fermi level) , which exists in a finite range between 
the metallic phase and the Mott-Hubbard-gapped phase 
present for half filling. 



Acknowledgments 

This work was supported in part by the Deutsche 
Forschungsgemeinschaft through SFB 608. 



1 For a review, see S. V. Kravchenko and M. P. Sarachik, 
Rep. Prog. Phys. 67, 1 (2004). 

2 P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

3 E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. 
V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 

4 B. Srinivasan, G. Benenti, and D. L. Shepelyansky, Phys. 
Rev. B 67, 205112 (2003). 

5 D. Heidarian and N. Trivedi, Phys. Rev. Lett. 93, 126401 
(2004). 

6 K. Kobayashi, B. Lee, and N. Trivedi, arXiv:0807.3372 
(unpublished) . 

7 P. B. Chakraborty, P. J. H. Denteneer, and R. T. Scalettar, 
Phys. Rev. B 75, 125117 (2007). 

8 Y. Song, R. Wortis, and W. A. Atkinson, Phys. Rev. B 77, 
054202 (2008). 

9 I. F. Herbut, Phys. Rev. B 63, 113102 (2001). 

10 D. Tanaskovic, V. Dobrosavljevic, E. Abrahams, and G. 
Kotliar, Phys. Rev. Lett. 91, 066603 (2003). 

11 M. C. O. Aguiar, V. Dobrosavljevic, E. Abrahams, and G. 
Kotliar, Physica B 403, 1417 (2008). 

12 K. Byczuk, W. Hofstetter, and D. Vollhardt, Phys. Rev. 
Lett. 94, 056404 (2005). 

13 E. Z. Kuchinskii, I. A. Nekrasov, and M. V. Sadovskii, 
JETP 106, 581 (2008). 

14 For recent references on that topic, see, e.g., D. Clement, 
P. Bouyer, A. Aspect, and L. Sanchez-Palencia, Phys. Rev. 
A 77, 033631 (2008); L. Fallani, C. Fort, and M. Inguscio, 
arXiv:0804.2888, Adv. At., Mol., and Opt. Phys. (to be 
published) . 

15 Such effect has also been observed in models with long- 
range Coulomb interactions. See Z. A. Nemeth and J.-L. 
Pichard, Eur. Phys. J. B 45, 111 (2005), and references 
therein. 

16 B. L. Altshuler and A. G. Aronov, in Electron- Electron 
Interactions in Disordered Systems, edited by A. L. Efros 
and M. Pollak (North Holland, Amsterdam, 1985), p. 1. 

17 P. Henseler, J. Kroha, and B. Shapiro, Phys. Rev. B 77, 
075101 (2008). 

18 D. J. Thouless, J. Phys. C 5, 77 (1972). 

19 D. Vollhardt and P. Wolfle, Phys. Rev. Lett. 45, 842 
(1980); Phys. Rev. B 22, 4666 (1980); Phys. Rev. Lett. 



48, 699 (1982). 

20 D. Vollhardt and P. Wolfle, in Electronic Phase Transi- 
tions, edited by W. Hanke and Y. V. Kopaev (Elsevier 
Science, Amsterdam, 1992), p. 1. 

21 T. Kopp, J. Phys. C 17, 1897 (1984); J. Phys. C 17, 1919 
(1984). 

22 A recent DMRG study of the ID Anderson-Hubbard model 
supports this assumption. See M. Okumura, S. Yamada, N. 
Taniguchi, and M. Machida, Phys. Rev. Lett. 101, 016407 
(2008). 

23 Note that in d = 3 £ diverges at the Anderson transition, 
and in d = 2 it can become large even for A/t ~ O(10) 3> 
1, as seen in Sec. IV. 

24 E. Medina and M. Kardar, Phys. Rev. B 46, 9984 (1992). 

25 J. Kroha, T. Kopp, and P. Wolfle, Phys. Rev. B 41, 888 
(1990). 

26 J. Kroha, Physica A 167, 231 (1990). 

27 B. Bulka, M. Schreiber, and B. Kramer, Z. Phys. B 66, 21 
(1987). 

28 Note the relative factor of 2 in the definitions of the decay 
length of the wave function and of the density correlations, 

( 6 = 2£- 

29 F. Yonezawa and K. Morigaki, Prog. Theor. Phys. 53, 1 
(1973). 

30 J. Hubbard, Proc. Soc. London, Ser. A 281, 401 (1964). 

31 B. Velicky, S. Kirkpatrick, and H. Ehrenreich, Phys. Rev. 
175, 747 (1968). 

32 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

33 Note the factor of 1/2 in the definitions of the disorder 
parameter in Ref. and the present study, V = A/2. 

34 P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

35 For a similar discussion in the case of a long-ranged inter- 
action see, e.g., G. Fleury and X. Waintal, arXiv:0807.3433 
(unpublished) . 

36 This value is somewhat smaller than the one obtained in 
Ref.l26l. because there the lattice periodicity of the particle- 
hole propagator in the center-of-mass momentum Q was 
not considered. 



